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We calculate within the Hard Thermal Loop expansion the produc- 
tion rate of soft photons and lepton pairs by a hot quark-gluon plasma 
at thermal equilibrium, up to the 2-loop order. Strong collinear di- 
vergences appear to mix the orders of the perturbative expansion so 
0^ ' that the 2-loop order is in fact dominant compared to the 1-loop order. 

More precisely, angular integrals that would otherwise be of order 1 
O |' are found to behave as powers of 1/g when the produced photon is 

JL . massless, therefore breaking the standard HTL power counting rules. 

S-H. 
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1 Introduction 

h ; 

The production rate of photons or dileptons is of some relevance to the phe- 
nomenology of a quark-gluon plasma. Indeed, these electromagnetic probes 
are interacting very slightly through electromagnetic forces. Therefore, it is 
expected that they will escape the plasma without re-interacting after their 
creation. As a consequence, they are rather clean probes of the state of the 
system at the time it produced them. 

Among the theoretical tools to calculate such observables, one can quote 
the semi-classical methods based on an eikonal approximation [1 — 4] , and 
the tools of finite temperature field theory [5 — 10]. In this report, I focus 
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mainly on the field theoretical method. I assume a plasma in thermal equi- 
librium and sufficiently hot so that the temperature is much larger than the 
masses of the particles. From the theoretical point of view, studying these 
production rates is a good test of our present understanding of problems 
such as infrared and collinear singularities in high temperature QCD. 

The starting point is to connect the photon production rate to the imag- 
inary part of the photon polarization tensor, via the following two relations 
[5,11]: 

Basically, these two formulae differ only by the allowed phase space, by the 
coupling constant involved in the decay of a virtual photon into a lepton 
pair and by the propagator of such a heavy photon 1 . 



2 Infrared problems and hard thermal loops 

2.1 Origin of the problem 

Before going on with some details of the calculation we performed, it is 
worth recalling some of the issues related to the concept of hard thermal 
loops (HTL) [12], and its relevance to the problem of infrared divergences. 
Indeed, although HTLs have been introduced to solve the problem of the 
apparent non gauge invariance of the gluon damping rate [13], they proved 
to be quite efficient in improving the infrared behavior of thermal gauge 
theories, which is a priori more severe than at T = due to the singular 
behavior of Bose-Einstein's factors at small energies: 

n B (/o) = eW r_ 1 «^»l in o «T. (3) 

Since in the Feynman rules, these factors are always accompanied by the 
Dirac's distribution 5(1%— l 2 — m 2 ), the Bose-Einstein's weights are a problem 
only if there are massless bosons in the theory. In QCD, one has gluons... 
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Figure 1: Hard and soft 
scales. 



2.2 Semantics: hard and soft scales 

It is useful in thermal field theories to distinguish between two typical energy 
scales: particles of energy of order T are said to be hard, and particles of 
energy of order gT are said to be soft. Roughly speaking, the hard scale is 
the typical energy of partons in the plasma, and the soft scale is the typical 
energy of quanta exchanged during parton interactions. 

2.3 Summation of hard thermal loops 

The foundations of the HTL concept lie in the fact that some one-loop dia- 
grams have a thermal contribution which is dominated by the hard region of 
phase space, and moreover this thermal contribution is as large as the bare 
corresponding function when all its external legs carry soft momenta. For 
instance, this is what happens for the one-loop photon self-energy in QED. 
Then, in order to take these large contributions into account, one defines 



an effective field theory [14, 15] by the means of effective propagators and 
vertices. The effective propagators are obtained by a Dyson summation of 
the HTL contribution to the corresponding self energy, whereas the effective 
vertices are defined to be the sum of their bare counterpart plus the HTL 
contribution to the corresponding vertex function. In order to avoid dou- 
ble counting of thermal corrections, one should include counter-terms whose 
purpose is to subtract at higher orders what had been added by this sum- 
mation. In that way, the overall Lagrangian function remains unchanged, 
so that the use of this effective theory just amounts to a re-ordering of the 
perturbative expansion. 

lr The above two formulae are true to all orders in the strong coupling constant a s , but 
only to the first non vanishing order in the QED coupling a, due to the fact that they 
neglect potential re-interactions of the produced photon on its way out of the plasma. 



Figure 2: Hard thermal loop 
of the photon self energy. 
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2.4 Basic properties 

The HTLs have a few nice properties that make the effective theories ob- 
tained by they summation rather attractive: 

• The HTLs are gauge invariant. As a consequence, the effective theory one 
obtains by summing all the HTLs is also gauge invariant. It is important to 
note the fact that the effective vertices are essential here to achieve the gauge 
invariance of the effective Lagrangian. In fact, starting by the summation 
of 2-point HTLs in order to obtain a better behaved effective propagator, 
one can obtain all the n— point (n > 3) HTLs by the requirement that the 
effective theory should be gauge invariant [14]. 

• The common physical interpretation of the HTLs is that they reflect the 
long distance behavior of the gauge interaction better than the bare theory 
did. In particular, they account for the Debye screening phenomenon, which 
is translated in field theory by the fact that the gluon acquires an effective 
mass m D ~ gT, thus leading to a screened interaction whose range is of 
order m~~ 1 . 



Figure 3: QED Debye screen 
ing. 

• The n— point (n > 3) HTLs are vanishing if one takes the trace over the 
Lorentz indices of two external gauge bosons. 

• Another property which makes the HTLs simple is the fact that the 1-loop 
diagrams that give the HTLs are calculated with massless particles running 
in the loop. Besides the technical simplification this fact provides, it is also 
at the origin of collinear singularities encountered when one is using HTLs 
with light-like external legs. 

2.5 Further problems 

Despite the great improvement the HTLs provide for perturbative calcula- 
tions in thermal gauge theories, some nasty problems remain. 

• The first one is the problem of collinear singularities evocated above. Due 
to the fact that the HTLs are calculated with massless propagators running 
in the hard loop, these functions are divergent when at least one of the 
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external legs is on-shell. This is due to an angular integral like 2 : 



+1 

dCl„ ( d cos 9 



77 « u„ ITTTn = 00 > ( 4 ) 



PK J k{l — cosQ) 
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where P = (l,p), and i^T is a light-like external momentum. A possible 
solution to this problem has been proposed recently [16], which amounts to 
keep an asymptotic thermal mass for the particle running in the hard loop. 
Obviously, the previous divergent integral becomes now: 



+i 

f d£l v f dcosO . , 

/ - °c / ; 77 < oo , (5) 

J P ■ K J Uk ~ k cos e 

-1 

where u>k = yj{k 2 + m?). Moreover, it has been shown that these improved 
hard thermal loops still generate a gauge invariant effective theory. 
• The other problem which is not solved by the summation of HTLs is due 
to the fact that static transverse gauge bosons don't get a thermal mass. 
This property is known to be true to all orders in QED, in which case it is 
related to the fact that static magnetic fields are not screened in a plasma. In 
QCD, there is still a hope that self-interactions of gluons can generate such 
a mass for the transverse modes, but it seems that this mass will be beyond 
the abilities of perturbative methods, and at most of order g 2 T contrary 
to usual thermal masses that are of order gT. The nullity of this magnetic 
mass is at the origin of further infrared divergences, like those encountered 
in the perturbative calculation of the damping rate of a fast fermion. 



3 7 production rate at 1 loop 
3.1 1-loop diagrams and results 

These 1-loop diagrams have been evaluated by several groups [5 — 8], leading 
to the following results: 

Imn%| [ai] (g ,q)~eV— , (6) 

ImlP H[6l] ((fo,g) = 0, Imn%| H (g o ,qf)=0. (7) 

2 This example displays only a logarithmic divergence. Nevertheless, when two external 
momenta are on-shell and collinear, one obtains stronger collinear singularities. 
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Figure 4: Contributions to the 
photon production rate at 1-loop in 
the effective theory. The first one is 
accompanied by the corresponding 
physical amplitudes. Black dots de- 
note effective propagators and ver- 
tices. 
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The result given for the diagram [ai] is to be preceded by a numerical 
coefficient which is finite if the emitted photon is massive (Q 2 > 0), but 
that diverges logarithmically when the photon mass goes to zero. This is 
due to the fact that the hard loop contained in the effective ^qq vertex is 
plagued by a collinear divergence when one of its external legs is light-like. If 
one applies the prescription presented in the previous section, this coefficient 
becomes a number of order ln(l/^). 



3.2 Comments 

A few remarks are worth saying concerning these results: 

• The only non-vanishing contribution at I-loop in the effective theory based 
on the summation of the HTLs comes from the diagram [ai] and corresponds 
to the interference between the two amplitudes represented in the figure 4. 
If, instead of using thermal field theory, one had just tried to guess what 
are the dominant processes contributing to photon production by a hot 
plasma, it is unlikely that these processes would have been thought as being 
dominant. 

• The exact nullity of the diagrams [6i] and [ci] is a consequence of the 
remark made earlier concerning the trace over Lorentz indices of n— point 
HTLs. 

• In the case of the diagram [ax], taking the trace ■ ■ doesn't make the 
result vanish, but nevertheless is at the origin of a suppression by a power 
of g. Therefore, if this suppression is specific to this topology, it may be 
that some 2-loop diagrams are of the same order (in thermal field theories 
the expansion parameter is g instead of g 2 ). 

This fact, added to the oddity of the dominant physical process at 1-loop, 
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should incite one to evaluate the 2-loop contributions. 



4 7 production rate at 2 loops 
4.1 2-loop diagrams 

In the list of 2-loop diagrams contributing to the photon production rate, 
we skipped right from the beginning those having a tadpole topology since 
taking the • • • M M trace will make them exactly cancel. Moreover, we skipped 
topologies containing effective vertices like "fggg or 'jgqq. This will be justi- 
fied later by the fact that the topologies [02] and [62] are important because 
of quasi superposed collinear singularities (leading to an enhancement by 
powers of l/g), whereas the diagrams we skipped can have at most powers 
of ln(l/p). For the remaining diagrams, it is a priori essential to keep the 
topologies with counter-terms in order to avoid double counting of thermal 
corrections already included at the 1-loop level. Nevertheless, it is obvious 
that the diagrams with counter-terms are of the order of magnitude of [ax] 3 
(the counter-term is part of the effective vertex contained in [ai]). 



Figure 5: Some 
2-loop contributions 
to photon production. 
The cross denote HTL 
counter-terms. 




4.2 Physical processes 

Looking more carefully at the kinematics of these 2-loop diagrams, we can 
see that the dominant contribution comes from the region of phase space 
where the quark running in the loop is hard while the exchanged gluon 
is soft. As a consequence, we can approximate the effective vertices by 
their bare counterpart, and the effective quark propagators by their hard 

3 If the effective theory works properly, the diagrams [<Z2] and [62] should be partly 
compensated by these counter-term diagrams, leaving us with a sub-dominant correction. 
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approximation: 



S(P) 



■ Vol 



Pi 



u\(p) + ie 



(8) 



with uj+(p) = \/(p 2 + M 2 ), where M 2 = g 2 C F T 2 /A. Therefore, we can 
replace the diagrams [a2] and [b2] by the simpler version of the figure 6, 
where we used a slightly abusive notation for the quark propagators we 
represented as bare ones whereas we are retaining an asymptotic thermal 
mass in them. On the same figure, we also represented the cuts relevant to 



Figure 6: Simplified 
version of [a2] and [b2]. 





the calculation of the imaginary part of the polarization tensor. The cuts 
not represented always imply that the quark be soft, which would drastically 
reduce the phase space available. A priori, by cutting an effective gluon 
propagator, we can pick the pole part (L 2 > 0) or the Landau damping part 
[1? < 0) of the gluon spectral function. The two contributions correspond 
respectively to photon production by Compton effect or by bremsstrahlung, 
which looks more intuitive than the processes involved in [al]. 




Figure 7: Physical processes involved in [b2]. 



4.3 Collinear behavior 

If both the external photon and the internal quark are massless, the diagrams 
[a2] and [b2] exhibits a collinear singularity when the 3-momentum of the 
quark is collinear to that of the external photon. Moreover, contrary to the 
case of the diagram [al] where such singularities are at most logarithmic, 
we have here two denominators (corresponding to the propagators which 
are not cut) that vanish almost simultaneously. This is obvious in the case 
of the diagram [a2] for which we have an exact double pole due to the 
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insertion of the self-energy correction. In the case of the diagram [b2], both 
denominators are distinct and so are the two poles. At Q 2 = 0, R 2 is 
vanishing when p is parallel to q, and (P + L) 2 is vanishing when r + 1 is 
parallel to q. Nevertheless, since p and r are hard whereas q and 1 are soft, 
these two conditions of collinearity are satisfied almost at the same time. It 
means that the two poles are very close, so that their association behaves 
almost like a double pole. 

Therefore, in both [a2] and [b2], we have denominators that lead to a lin- 
ear collinear divergence. A careful calculation of the angular integral while 
keeping Q 2 > and the quark thermal mass M F indicates that these diver- 
gences are regularized by the combination M 2 S = M 2 + Q 2 r 2 /q 2 divided by 
T 2 . If we limit ourselves to photons of very small virtuality (Q 2 < then 
this dimensionless regulator is much smaller than 1, and the dimensionless 
angular integral is much larger than 1, which would be its order of magni- 
tude in the absence of collinear problems. In fact, at Q 2 = 0, the only scale 
in the angular integral is M 2 /T 2 ~ g 2 , so that the enhancement factor of 
this integral is l/g 2 . The fact that powers of 1/g can appear at Q 2 = in 
the perturbative expansion due to collinear singularities can make the order 
of magnitude of a diagram somewhat unpredictable by the power counting 
rules associated usually with the HTL expansion (indeed, these rules always 
assume that the dimensionless angular integrals that one can factorize out 
of the HTLs are of order 1, and we have shown by this example that these 
integrals can sometimes be as large as powers of l/g). 

The above results concerning the order of magnitude of the angular inte- 
grals indicate that the result is completely dominated by the collinear sector 
of the phase space. As a consequence, one can use a collinear approximation 
to obtain a simplified expression for the Dirac's algebra associated with the 
quark loop. The results one obtains that way are: 

Tr [a2] « l(yQ 2 P a Pf3 , Tr [b2] « 16(Q 2 - L 2 )P a P p , (9) 

where a and (3 are the Lorentz indices of the gluon. Then, limiting ourselves 
to very small Q 2 , we can neglect the contribution of the diagram [a2] . 

Moreover, since we have the constraint 5(P 2 — M 2 ), the denominator 
{P+L) 2 —M 2 can approach zero only if L 2 is negative. As a consequence, the 
portion of phase space that contains the collinear singularities also verifies 
L 2 < 0, which indicates that the bremsstrahlung production of photons 
dominates with respect to the Compton effect. 
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4.4 Results 

Therefore, the result we obtain for the imaginary part of the polarization 
tensor of the photon can take the form [9, 10]: 

ImII%(<Z ,q) « - ° 2 F E J t,l — , (10) 

where the quantities J T L are pure numbers quantifying the respective con- 
tributions of transverse and longitudinal gluons exchange. These coefficients 
are given by an integral which is evaluated numerically, as a function of the 
ratio Q 2 jq 2 and plotted on the figure 8. We can see clearly an enhancement 
in the region of small Q 2 . Indeed, since the regularization of collinear diver- 
gences is realized by the mass M 2 S , increasing with Q 2 , the result is larger 
at small Q 2 . The same quantities are plotted also with a logarithmic scale 



Figure 8: Transverse and longitu- 
dinal contributions as a function of 
Q 2 /qo, for 3 colors and 3 light fla- 
vors. The value taken for the cou- 
pling constant is g = 0.44. 



that enables one to see in more detail the region of small Q 2 . In particular, 
we see that they converge to finite values when Q 2 goes to 0. Moreover, we 
see that the enhancement effect (quantified by the ratio of the values taken 
at Q 2 = and at Q 2 ~ q 2 ) is larger when g is small. 



Figure 9: Effect of the value of 
the coupling constant. Solid lines: 
g — 0.1, dashed lines: g = 0.44, 
dotted lines: g = 2. We used 3 
colors and 3 light flavors. 



Let us now list a few important features of the result we obtained: 

• The diagram [al], only non vanishing contribution at the 1-loop order, 

is negligeable compared to [b2], for small enough Q 2 (Q 2 <C q 2 <C T 2 ). 
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Moreover, the latter is much larger than the diagrams with counter-terms 
of the figure 5. This fact is due to strong collinear divergences. Therefore, 
there is a possibility that such collinear singularities may invalidate the 
perturbative expansion since powers of the inverse coupling constant emerge 
in the process. 

• The result we obtained is free of any infrared divergence and is totally 
insensitive to the magnetic mass scale as long as the magnetic mass is neg- 
ligeable in front of the thermal mass M F . If this condition is satisfied, then 
the regularization of potential IR singularities in the transverse gluon con- 
tribution is done by the quark mass M F . There is a natural explanation for 
this a priori surprising property related to the fact that the delta function 
constraints become incompatible in the limit L — > when M F > 0. 

• Our result can be put into a form that looks like the expressions obtained 
by semi-classical methods: 



The interpretation of the various factors entering the above formula is rather 
easy: starting from the beginning, we find the photon phase space, the phase 
space of two incoming quarks, the phase space of the same outgoing quarks, 
the amplitude \ Ai\ 2 of the scattering process between the two quarks, associ- 
ated to a delta constraint ensuring the overall momentum conservation, and 
finally a factor which is the square of the electromagnetic current coupling 
the photon to the quark line. 

Therefore, we see that thermal field theory gives the same kind of expres- 
sion as the one obtained by semi-classical methods [1 — 4] when considering 
the photon emission induced by a single scattering. In fact, the thermal 
field theory result is more precise in the sense that it has not neglected the 
contribution of transverse gluons, which is as important as the longitudinal 
one as can be seen from the curves plotted above. 



dN 




(11) 
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5 Conclusions and perspectives 



The study of the photon production rate by a hot quark-gluon plasma by 
using the effective theory based on the summation of the hard thermal loops 
has shown that serious problems may arise from strong collinear singular- 
ities. More precisely, angular integrals that are usually assumed to be of 
order unity in the HTL framework are found to behave as powers oil/g. It 
is not known by now if such singularities or even stronger ones may emerge 
from higher order contributions, enabling them to remain at the same level. 

Prom a more phenomenological point of view, we arrived at the con- 
clusion that real or slightly virtual soft photons are mostly produced by 
bremsstrahlung in a hot plasma. 

An important issue is to determine whether higher order contributions, 
corresponding to photon emission induced by multiple scatterings, are im- 
portant or not. The importance of this study is related to the fact that the 
semi-classical methods predict that multiple scatterings are important and 
can modify the photon spectrum in the region of small energy, a phenomenon 
known as the Landau-Pomeranchuk-Migdal effect. 
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